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We compute the colour fields of SU(3) lattice QCD created by static pentaquark systems, in a 
24 3 x 48 lattice at /3 = 6.2 corresponding to a lattice spacing a — 0.07261(85) fm. We find that 
the pentaquark colour fields are well described by a multi-Y-type shaped flux tube. The flux tube 
junction points are compatible with Fermat-Steiner points minimizing the total flux tube length. 
We also compare the pentaquark flux tube profile with diquark-diantiquark central flux tube profile 
in the tetraquark and the quark-antiquark fundamental flux tube profile in the meson, and they 
match, thus showing that the pentaquark flux tubes are composed of fundamental flux tubes. 



INTRODUCTION 

Here we study the colour field flux tubes produced 
by static pentaquarks in SU(3) lattice QCD. Unlike the 
colour fields of simpler few-body systems, say mesons, 
baryons and hybrids, [1-4], the pentaquark fields have 
not been previously studied in lattice QCD. This study 
is relevant both for the solution of theoretical problems 
and for the development of phenomenological models of 
QCD. 

Quark confinement remains one of the main open theo- 
retical problems of particle physics. In lattice QCD, flux 
tubes composed of colour-electric and colour-magnetic 
fields have been observed and this constitutes a very im- 
portant clue for the understanding of quark confinement. 
Since the onset of QCD with its asymptotic freedom and 
infrared slavery, it is well known that confinement is due 
to the gluon fields and suppressed by the quark fields. It 
is thus important to measure the different possible flux 
tubes of pure gauge lattice QCD, to provide data for any 
theoretical attempt to solve the QCD confinement prob- 
lem. 

Moreover, in what concerns phenomenology, the study 
of the colour fields in a pentaquark is important to 
discriminate between different multi-quark Hamiltonian 
models, quark models with two-body interactions only [5] 
as in the original quark model, from the string flip-flop 
model with a multi-body potential [6] . In the string flip- 
flop model, the colour charges are connected by strings 
disposed geometrically in order to minimize the total 
string length. The strings constitute the limit of very 
thin elementary flux tubes. An elementary or fundamen- 
tal flux tube is the flux tube connecting the quark and 
antiquark of a meson, where the quark is in the triplet or 
fundamental representation of QCD. For instance in the 
two quark two antiquark system, depending of the posi- 
tion of these colour charges, the minimal string may be 
a two-meson string, or a tetraquark string, shaped like a 
double-Y flux tube, as in Fig. 1, composed of five linear 
fundamental flux tubes meeting in two Fermat-Steiner 
points [7-9]. A Fermat, or Steiner, point is defined as a 




Figure 1: (Colour online.) In the string flip-flop model, 
thin elementary flux tubes similar to strings connect the 
colour charges in order to minimize the total length of 

the strings. Whenever geometrically possible, three 
elementary flux tubes meet in a Fermat-Steiner point at 
an angle of a — 120°. Here we depict planar examples 
of a tetraquark and a pentaquark flux tubes. 



junction minimizing the total length of strings, where lin- 
ear individual strings join at 120° angles. When the posi- 
tions of the colour charges change, the potential may thus 
flip from one four-body potential to a pair of two-body 
potentials and flop back again. Notice the flip flop poten- 
tial, compatible with the confining component of the flux 
tubes explored here, lead to tetraquark boundstates, be- 
low the strong decay threshold to pairs of mesons [7, 10- 
12]. Recent investigations found that, even above the 
strong decay threshold, the presence of a centrifugal bar- 
rier in high angular momentum multiquarks may increase 
the stability of the system [6, 13]. The multiquark hamil- 
tonians are important to understand not only the elusive 
multiquark hadrons, but also high density QCD where 
many quarks may overlap. 

Experimentally, multiquark exotic hadrons have been 
searched for many years because as soon as the quarks 
were proposed in the sixties to classify the meson and 
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Figure 2: Pentaquark Wilson loop as denned by 
Okiharu et al. [15, 16]. Here we extend this Wilson loop 
with different paths Li and Ri for the quarks and Mi 
for the antiquark. 



baryon resonances, and the quark model was proposed in 
the seventies [14], it became clear that systems with more 
than three quarks could also possibly exist. One of the 
main problems of hadronic physics is thus to determine 
whether multiquark resonances exist or not, and weather 
the possible multiquark resonances are narrow or wide. 

The simplest multiquark system is the tetraquark, and 
it was already proposed by Jaffe in the seventies [17] as 
a bound state formed by two quarks and two antiquarks. 
Presently some observed resonances are tetraquark can- 
didates. Very recently the Belle Collaboration made the 
tantalizing observation [18], in five different T(5S) decay 
channels of two new charged bottomonium resonances 
Z b with masses of 10610MeV/c 2 and 10650MeV/c 2 and 
narrow widths of the order or 15 MeV. Since all stan- 
dard bottomonia are neutrally charged, these two new 
resonances have a flavour only compatible with bbud 
tetraquarks. In 2003, the X(3872) observed by the Bell 
Collaboration [19, 20] was suggested as a tetraquark can- 
didate by Maiani et al [21]. In 2004, the D sJ (2632) state 
seen in Fermilab's SELEX [22, 23] was suggested as a pos- 
sible tetraquark candidate. In 2009, Fermilab announced 
the discovery of Y(4140), which may also be a tetraquark 
[24] . There are as well indications that the Y(4660) could 
be a tetraquark state [25]. The T(5S) bottomonium has 
also been recently suggested to be a tetraquark resonance 
[26]. However a better understanding of tetraquarks is 
necessary to confirm or disprove the X, Y, Z and possi- 
bly also other light resonances candidates as tetraquark 
states. 

The pentaquark is the next in the multiquark hadron 
series. Pentaquark hadrons were already proposed in the 
eighties by Manohar [27] and Chemtob [28], inspired by 
extensions of the Skyrme model. In the 2000s a burst of 
interest was sparkled by a discovery claim of the 9 pen- 
taquark by Nakano et al [29]. This led to many experi- 
mental and lattice QCD studies of pentaquarks, together 
with hundreds of theoretical estimations of the O prop- 
erties. However the resonance 6 ended up by not being 
confirmed by the scientific community [30, 31]. The many 
hundreds of publications on the subject, with disparate 
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Figure 3: Projection in the spacial dimensions of the 
different Wilson loop geometries for the static 
pentaquark studied in this work. The solid dots 
correspond to the quarks positions and the open dots to 
the antiquarks. The solid lines correspond to the 
space-like Wilson paths. 



conclusions, show that the O pentaquark was beyond the 
scope of the scientific techniques utilized in the 2000s. 

The multiquark hadrons are thus very elusive systems, 
much harder to observe experimentally, to understand in 
models, and to simulate in lattice QCD than the con- 
ventional mesons and baryons. Nevertheless, inasmuch 
as the understanding of confinement, the existence/non- 
existence of multiquark hadrons remain an important 
problem in QCD, to be further explored in the future 
PANDA experiment at GSI. 

It is thus important to proceed with the well defined 
program of understanding the static potentials and flux 
tubes of multiquarks in quenched Lattice QCD. 

In the last years, the static tetraquark potential 
has been studied in Lattice QCD computations [32- 
34]. The authors concluded that when the quark-quark 
are well separated from the antiquark-antiquark, the 
tetraquark potential is consistent with One Gluon Ex- 
change Coulomb potentials plus a four-body confining 
potential, suggesting the formation of a double-Y flux 
tube, typical of the four-body potential of the string flip- 
flop model as in Fig. 1, composed of five linear funda- 
mental flux tubes meeting in two Fermat-Steiner points 
[7—9]. This flux tube geometry was confirmed by Lat- 
tice QCD studies of the flux tubes produced by a static 
tetraquark system [35, 36]. In what concerns the pen- 
taquark, static potentials have already been explored in 
a geometry with the antiquark situated in the centre of 
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Table I: Pentaquark geometries studied and number of 
lattice configurations used in this work. The geometry 
type is outlined in Fig. 3. The column Id corresponds 
to the numbering used in the text. 



Id 


Geometry Type 


di 


d 2 


d 3 


# Configs. 


(i) 


Fig. 3a 


8 


8 





551 


(") 


Fig. 3a 


8 


8 


6 


549 


(iii) 


Fig. 3b 


4 


4 


8 


544 


(iv) 


Fig. 3b 


6 


4 


8 


1121 



Table II: Fermat-Steiner points for the pentaquark 
geometries studied. The geometry type is outlined in 
Fig. 3 and Table I. 

Fermat-Steiner Points 
Id n m r m 

(i) (0,-1.691,0) (0,1.691,0) (0,0,0) 

(ii) (3.897,-3.830,0) (3.897,3.830,0) (6,0,0) 

(iii) (-2.309,-4,0) (-2.309,4,0) (0,0,0) 

(iv) (-2.309,-6,0) (-2.309,6,0) (1.155,0,0) 



the four quarks [15, 32], also consistent with a string flip- 
flop model, in this case with only six fundamental flux 
tubes and two Fermat-Steiner points. 

Here we proceed with the flux tube research program, 
studying the flux tubes of static pentaquarks in pure 
gauge SU(3) lattice QCD. In Section II we detail the 
framework we set to measure the flux tubes. We also ex- 
tend the geometries explored in the static potential stud- 
ies. In Section III we expose our results and conclude. 



SIMULATING THE PENTAQUARK FLUX 
TUBES IN LATTICE QCD 

The static potential for the pentaquark was already 
studied in the lattice QCD by [32] and [15, 16] utilizing 
generalized Wilson loops. Here we use similar Wilson 
loops to place a static system of four quarks and one an- 
tiquark in the lattice, in four different geometries. More- 
over we measure the colour-electric and colour-magnetic 
fields produced by the static charges. 

The Wilson loop operator for the pentaquark system 
is defined in a gauge-invariant way, as illustrated in Fig. 
2, by 

W 5Q = A^V^'M"' (R 3 R 12 R 4 ) jj ' (L 3 L 12 L 4 ) kk ' , 



With this definition, the chromofields are given by 



where 



R i/i = l e ijk e i'j'k' R jj' R kk' ^ 

h{\ = ^e ijk (?' j ' k ' L{ j ' Lf' . 



(1) 
(2) 



The projection in the spacial dimensions of our four dif- 
ferent Wilson loop geometries for the static pentaquark 
is illustrated in Fig. 3. The distances of the geometries 
are quantified in Table I. The labelling of the geometries 
together with the number of lattice configurations used 
in this work are also shown in Table I. 

We compute the colour electric and the colour mag- 
netic fields, by using the correlators of the plaquettes 
P^ v and the Wilson loop W$q. We define the plaquettes 
as = 1 - iTr[[/ M (s)^(s + M )[/J(s + v)U}(b)]. 



(E 2 ) = (P 0i ) 



(W 5Q Ppj) 
(W 5Q ) 



jk) 



(W 5Q ) 



(Pjk) 



(3) 
(4) 



with the indices j and k complementing index i. The 
lagrangian and energy densities are given by C = \{E 2 — 
B 2 ) and % = \{E 2 + B 2 ). 

To compute the static field expectation value, we plot 
the expectation value (£f(r)) or (_B 2 (r)) as a function 
of the temporal extension T of the Wilson loop. At suf- 
ficiently large T, the groundstate corresponding to the 
studied quantum numbers dominates, and the expecta- 
tion value tends to a horizontal plateau. In order to im- 
prove the signal to noise ratio of the Wilson loop, we 
use 50 iterations of APE Smearing with w — 0.2 (as in 
[3, 35]) in the spatial directions and one iteration of hy- 
percubic blocking (HYP) in the temporal direction, [37], 
with ct\ = 0.75, a 2 = 0.6 and a 3 — 0.3. Note that these 
two procedures are only applied to the Wilson Loop, not 
to the plaquette. 

To check if the pentaquark flux tube produces a 
clear signal, we study the x 2 /dof of our pentaquark T 
plateaux. But, surprisingly, event at some of the dis- 
tances illustrated in Fig. 4, where the string flip-flop 
potential would favour the meson-baryon flux tube, with 
a lower energy than the pentaquark flux tube, we find 
T plateaux with a good x 2 /dof. This shows that the 
mixing between the pentaquark flux tube and the meson- 
baryon flux tube is small, and it is possible to study clear 
pentaquark flux tubes even at relatively large diquark 
distances. 

To compute the fields, we fit the horizontal plateaux 
obtained for each point r determined by the plaquette 
position, but we consider z — for simplicity. We finally 
compute the error bars of the fields with the jackknife 
method. 

We compute the Fermat-Steiner points with the it- 
erative method of Bicudo et al. [8]. We have five 
quarks (antiquarks) with the label i and three Fermat- 
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Figure 4: (Colour online.) Density plots of the chromoelectric and chromomagnetic fields and Lagrangian and 
energy densities for the geometries defined in Table I. The black dot points correspond to the Fermat-Steiner points. 

Table II. The results are presented in lattice spacing units. 



Steiner points with label a — I, II, III, Fermat-Steiner point, 



(Xi , Hi , %i ) ; 
(x a ,y a ,Z a ^ ) 

r la = \J {X a ~ Xi) 2 + (y a - Hi) 2 + (z a - Zi) 2 . (5) 



To minimize the total length of the strings. 

RESULTS AND CONCLUSION 

d = j»i i + r 2 i +r 3 u + r iU + rg m + rj m + r H ill , (6) 

We remark that the signal is clear only if the paths con- 
sidered in the Wilson loop overlap the flux tube. Thus we 
we only need to solve one non-linear vector equation per consider geometries for the Wilson loop where the paths 
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Figure 5: (Colour online) we show three-dimensional plots of the Lagrangian density for the geometries defined in 
Table I. The density enhancement is maximal in the location of the colour charges, and if the colour charges were 
close the Coulomb potential would be important. With our geometries the colour charges are separated, and the 
fundamental flux tubes connecting the charges and the Fermat-Steiner points are evident. The results are presented 

in lattice spacing units. 



are just some lattice spacings distant from the expected 
string position in the string flip-flop model. We consider 
the four different Wilson loop geometries, detailed in Fig. 
3 and in Table I. We only utilize planar geometries for 
the colour sources, in order to produce clearer pictures 
of the fields. The results for the colour field densities are 
presented only for the x y plane since the colour sources 
are in this plane and the results with z ^ are less in- 
teresting for this study. Then with colour field densities 
as a function of x and y we produce density plots and 
three-dimensional plots. 

To produce the results presented in this work , we 
utilize quenched configurations in a 24 3 x 48 lattice at 
(3 — 6.2. The number of configurations used is described 
in Table I. We present our results in lattice spacing units 
of a, with a = 0.07261(85) fm or a" 1 = 2718 ± 32 
MeV. We generate our configurations in NVIDIA CPUs 
of the FERMI series (480, 580 and Tesla C2070) with 
a SU(3) CUDA code upgraded from our SU(2) combi- 
nation of Cabibbo-Marinari pseudoheatbath and over- 
relaxation algorithm [38-40]. Our SU(3) updates involve 
three SU(2) subgroups, we work with 9 complex numbers, 



and we reunitarize the matrix. 

The results for the colour fields, the energy and la- 
grangian densities are shown in Figs. 4—5. The figures 
clearly exhibit multi-Y-type shaped flux tubes. We also 
plot the Fermat-Steiner points defined in Table II. The 
Fermat-Steiner points of geometries (i) and (ii) are of dif- 
ferent type from the Fermat-Steiner points of geometries 
(iii) and (iv), since in the first geometries angles of 120° 
between the fundamental strings are not possible and 
thus the central Fermat-Steiner point has merged with 
the antiquark source. Nevertheless, and although the 
flux tubes have a finite width and are not infinitely thin 
as is assumed in the string flip-flop models, and although 
the Coulomb component of the potential is certainly im- 
portant, we notice the junctions for the elementary flux 
tubes are clearly close to the computed Fermat-Steiner 
points. This validates the use of string flip-flop models 
for the quark confinement in constituent quark models. 

In Fig. 6, we compare the chromoelectric field profile 
for the pentaquark, tetraquark and the quark-antiquark 
system in the middle of the flux tube. The tetraquark 
and the quark-antiquark results were obtained by [35]. 
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Figure 6: (Colour online.) Profile of the chromoelectric 
field for the pentaquark, tetraquark and 
quark-antiquark systems. The pentaquark profile 
corresponds to the geometry profile outlined in Fig. 3 

and Table I along x = 4. The tetraquark and 
quark-antiquark results are from [35], in the middle of 
the flux tube. 



The three chromoelectric fields are identical up to the er- 
ror bars. This confirms that the pentaquark flux tube 
is composed of a set of fundamental flux tubes with 
Fermat-Steiner junctions, and again validates the string 
flip-flop models as models for the quark confinement in 
constituent quark models. 

Multiquark stability is a subtle theoretical problem, 
requiring the correct understanding and calibration of 
the quark interactions. Combining our pentaquark re- 
sults with the flux tube studies of mesons [41] , baryons 
[42], hybrids [3], glueballs [43], and tetraquarks [35, 36] 
we finally feel confident that the string flip-flop poten- 
tial, where fundamental strings with the minimal possi- 
ble length link the static colour sources, is the correct 
phenomenological model for the confinement of any sys- 
tem of static quarks, antiquarks and gluons. Whether 
the string flip-flop confining potential together with a 
correct short-range potential lead to multiquark narrow 
resonances or boundstates remains a difficult quantum 
mechanical problem, but very interesting to the confine- 
ment and quark model experts. 
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